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Spin tunneling in molecular magnets has attracted much 
attention, however theoretical considerations of this phe- 
nomenon up to now have not taken into account the many- 
spin nature of molecular magnets. We present, to our knowl- 
edge, the first successful attempt of a realistic calculation 
of tunneling splittings for Mni2 molecules, thus achieving a 
quantitatively accurate many-spin description of a real molec- 
ular magnet in the energy interval ranging from about 100 K 
down to 10~^^ K. Comparison with the results of the stan- 
dard single-spin model shows that many-spin effects affect 
the tunneling splittings considerably. The values of ground 
state splitting given by single-spin and many-spin models dif- 
fer from each other by a factor of five. 
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Progress in coordination chemistry lead recently to 
synthesis of a completely new class of magnetic mate- 
rials of nanometer size, molecular magnets, which has 
drawn the attention of physicists as well as chemists [Q . 
In particular, these materials have been proven to be very 
suitable for the study of mesoscopic quantum tunneling 
effects in magnetic materials. A number of impressive 
experimental results have been obtained recently, such 
as the thermally-induced tunneling ||^ , ground state - to 
- ground state tunneling and topological phase effects in 
spin tunneling ||^. 

Among others, the 

molecular magnet Mni20i2(CH3COO)i6(H20)4 (below 
referred to as Mni2) constitutes a subject of great in- 
terest. The effect of resonant magnetization tunneling 
has been first observed and studied in detailed experi- 
ments on Mni2 ||], and, at present, a substantial amount 
of reliable experimental data has been collected. How- 
ever, progress in understanding the properties of molec- 
ular magnets is greatly hampered by the lack of a com- 
prehensive theoretical description. Recent experiments 
show that the conventional single-spin description of 
magnetic molecules, including Mni2, does not accurately 
describe their properties and the constituent many-spin 
nature of these molecules should be taken into account. 
This can be particularly important for understanding 
tunneling properties; for which it has been demonstrated 
that the single-spin description of a many-spin sys- 
tem can give totally misleading results for tunneling split- 
tings, differing by three orders of magnitude from exact 
values. 

Therefore, the development of an adequate many-spin 
description of molecular magnets constitutes an impor- 
tant theoretical problem, which becomes especially com- 



plicated if the tunneling properties are to be studied. 
For example, the magnetic cluster of Mni2 molecule con- 
sists of eight ions Mn^+ with spin 2 and four Mn'*+ ions 
with spin 3/2 (see Fig. ^), coupled by anisotropic inter- 
actions, so that the Hilbert space of the corresponding 
spin Hamiltonian consists of 10^ levels. Furthermore, 
the tunneling splittings in Mni2 are very small, of or- 
der of 10^^^ Kelvin. Obviously, the brute-force direct 
calculation of tiny tunneling splittings, even for several 
low-lying states, is far beyond the capabilities of modern 
computers. Although several methods have been devel- 
oped to cope with this problem j^,^ , we are not aware of 
any calculations of tunneling splittings made for realistic 
models of complex molecular magnets, so that the pos- 
sibility of such a simulation has not been demonstrated 
yet. 

In this paper, we present, to our knowledge, the first 
successful attempt to calculate the tunneling splitting for 
a realistic model of Mni2 [Q, thus explicitly showing the 
possibility of calculating tunneling splittings in rather 
complex molecular magnets. As a result, we achieve a 
quantitatively accurate many-spin description of a real 
system in the energy interval ranging from about 100 K 
down to 10^1^ K, covering 14 orders of magnitude. 

The 8-spin model, which forms a basis for present cal- 
culations, has been described and discussed thoroughly 
in Ref. 0] . The central idea of the model is to use the nat- 
ural hierarchy of interactions present in Mni2. Namely, 
the antiferromagnetic exchange interactions Ji ~ 220 K 
(see Fig. ^ between ions Mn^+ and Mn''+ Q are sig- 
nificantly stronger than other isotropic exchange (J2, J3 
and J4) and anisotropic interactions, so at low tempera- 
tures they can be considered as absolutely rigid. In other 
words, the spins of the corresponding ions can be consid- 
ered as forming dimers with a total rigid spin 1/2. The 
validity of this model has been supported by megagauss- 
field experiments | pO[ : the states of dimers with the spin 
higher than 1/2 (excitations of dimers) come into play 
when the external magnetic field is about 400 T. The 
same conclusion can be drawn from the the dependence 
of the magnetic susceptibility of Mni2 versus tempera- 
ture (see, e.g. Ref. Q). 

After consideration of different possible interactions, 
the following spin Hamiltonian has been proposed in Ref. 
to describe Mni2 molecules: 
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Here, Sj and are the spin operators of large spins 
S = 2 and small dimer spins s = 1/2, correspond- 
ingly, where the subscript i denotes the index of the 
spin. Isotropic exchange between small spins and large 
spins is described by the parameter J', whereas J denotes 
the exchange of small spins with each other. The third 
term describes the single-ion uniaxial anisotropy of large 
spins. Finally, the fourth term describes the antisym- 
metric Dzyaloshinsky-Morya (DM) interactions in Mni2, 
and D*'-' is the Dzyaloshinsky-Morya vector describing 
DM-interaction between i-th small spin and j-th large 
spin. The molecules of Mni2 possess the fourfold rotary- 
reflection axis (the symmetry C4), which imposes restric- 
tions on the DM- vectors D*'^, so that Dzyaloshinsky- 
Morya interactions can be desribed by only three param- 
eters = Dy = Dl'^, and D, = D^^ . 

The model discussed above has been demonstrated to 
describe correctly the energy spectrum up to about 100 
K @,|ll[l. It has allowed explaination of a wide range of 
experimental data, such as the unexpected splitting of 
neutron scattering peaks [Q , results of EPR experiments 
and measurements of magnetic susceptibility versus 
temperature Q. This model has provided a quantita- 
tive description of the response of Mni2 molecules to the 
transversal magnetic field 1 11 1 (external field applied per- 
pendicular with respect to the easy axis of the molecule). 
The following set of the parameters has been determined 
for the spin Hamiltonian (0): 



J = 3.6 K, J' = 84 K, = 5.69 K 



(2) 



= 25.3 K, Dy 
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Having this model at hand, can we describe correctly the 
tunneling properties of Mni2? 

In the Hamiltonian (|^), only the fourth term, rep- 
resenting Dzyaloshinsky-Morya interactions, allows for 
tunneling. Indeed, the first two terms (describing 
isotropic exchange interactions) conserve both the total 
spin of the molecule S and its projection Sz (the z-axis 
is chosen to coincide with the C4 axis of the molecule). 
Thus, these terms can not lead to tunneling between the 
levels with different Sz ■ The third term, representing an 
easy-axis anisotropy, also conserves this quantity. There- 
fore, even though the levels \Sz = +M) and \Sz = ~M) 
are degenerate, in the absense of the Dzyaloshinsky- 
Morya term no tunneling between them can appear. But 
the DM-interaction mixes levels with different 5z, thus 
giving rise to tunneling of the molecule's spin. Below, for 
simplicity, we will denote the energy levels by the value of 
Sz- Although it is not an exact quantum number in the 
system under consideration, we can formally consider the 
DM-interaction as a perturbation, and use perturbation 
theory terminology. 



The first question to pose concerns the precision of 
the level splitting calculation. Parameters of the Hamil- 
tonian are determined with some finite precision, and a 
small error (say, of the order of several Kelvin) affects 
the level energy by an amount of order of Kelvin, which 
is much larger than the very small value of tunneling 
splitting (of order of 10~^^ K). Does it deprive the calcu- 
lational results of all meaning? To answer this question, 
we note that the levels \Sz = +M) and \Sz ~ —M) are 
degenerate due to exact symmetry properties of the spin 
Hamiltonian, and, in the absense of the DM-term, would 
be degenerate at any values of parameters. Therefore, the 
tunneling splittings IS.E^m.~m are governed only by the 
strength of the interaction which breaks the symmetry, 
i.e. DM-interaction. Therefore, if the parameters of the 
Hamiltonian are determined with reasonably small rel- 
ative error, and if the simulation is done with sufficient 
accuracy, then the relative error of the level splittings will 
also be small. The results of our calculations confirm this 
conclusion. 

In the calculation scheme we employed 0, the most 
significant error comes from the fact that only a fraction 
of all the levels produced by the spin Hamiltonian (|^) are 
used in the calculations, and the influence of levels having 
higher energy is neglected. Therefore, to achieve reason- 
able accuracy, a sufficiently large number of levels should 
be taken into account. We studied the dependence of 
the resulting tunneling splittings IS.E^m,-m for different 
pairs of degenerate levels \Sz = -\-M) and \Sz = —M) on 
the number of lowest levels actually used in calculations. 
The results are shown on Fig. |^. It can be seen that the 
reasonable accuracy can be achieved by accounting for 
about the 700 lowest levels. 

Final results for the tunneling splittings /S.E-\-m.-m 
are presented in the Table, where the values obtained 
in many-spin calculations are compared with those given 
by the single-spin model. It is important to mention 
that the tunneling splittings should be zero for the lev- 
els with odd values of M (i.e., \Sz = ±9), \Sz = ±7) 
etc.) because the fourfold symmetry of the molecule im- 
poses certain restrictions on the symmetry of the spin 
Hamiltonian and makes some matrix elements vanish. 
In the single-spin model of Mni2 this property of the 
spin Hamiltonian is introduced explicitly, by retaining 
only those operators which possess the required four- 
fold symmetry. In the many-spin simulations, we obtain 
the same result independently: the energies of the levels 
\Sz = ±M) for odd M are the same with the accuracy 
of order of 10^^^ K, i.e. of order of computational error. 
This value is much less than the smallest of the splittings 
A£;+io,-io = 2.03-10-12 K. 

Analyzing the results presented in the Table, we see 
that the single-spin and many-spin models give rather 
close results for higher energy levels M — ±2 and M — 
±4. It has been shown that these values of tunneling 
splittings allow to describe the experimental data avail- 
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able now with a good precision. On the other hand, for 
the ground-state levels M = ±10 the single-spin model 
predicts the value of tunneling splitting five times larger 
than the many-spin model. The possibility of such a dif- 
ference has been predicted before and, to our knowl- 
edge, Mni2 is the first example of a real system where 
this difference has been found. Unfortunately, reliable 
experimental data concerning ground state - to - ground 
state tunneling in Mni2 are absent, so the single-spin 
and many-spin models can not be distinguished on the 
basis of experimental results. Nevertheless, many-spin 
consideration can be important for a correct quantitative 
description of the ground state - to - ground state tunnel- 
ing in other molecular magnets, such as Feg. Many-spin 
effects may be one possible explanation for the disagree- 
ment found in Rcf. between the experimental results 
and predictions of the single-spin model. 

Summarizing, we evaluated the tunneling splittings in 
Mni2 on the basis of the realistic many-spin model. We 
demonstrated that even tiny energy splittings, of order of 
10^^'^ K, can be calculated with reasonable precision for a 
rather complex molecular magnet. The results obtained 
have been compared with predictions of the single-spin 
model. We found that both models give close values for 
splittings of higher levels. Thus, the models can not be 
distinguished on the basis of experimental results avail- 
able now, which provide the information only about up- 
per levels splittings (M = ±4 and M = ±2). How- 
ever, the ground state splittings calculated using these 
two models differ by a factor of five. This difference may 
be important also for other molecular magnets, such as 
Fes, where the ground state splitting given by the single- 
spin model is three times less than the value obtained in 
experiment. 
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FIG. 1. Schematic plot of the Mni2 cluster. Small black 
circles represent Mn*^ ions, large white circles — Mn^^ ions. 
Different types of lines connecting the ions (solid, dashed, dot- 
ted and dash-dotted) correspond to different types of isotropic 
Heisenberg exchange interactions (Ji, J2, J3 and J4). 

FIG. 2. Dependence of the logarithm of tunneling splittings 
AE+M,-M versus the number of levels taken into account in 
the many-spin calculations. The results for M — 10, 8, 6, 4, 
and 2 are presented. Tunneling splittings for the levels with 
odd M are zero, because of the symmetry properties of the 
spin Hamiltonian. 



TABLE 


I. Calculated values 


of tunneling splittings 


AE+M,-M 


; comparison between 


the single-spin and the 


many- spin 


models. 




M 


Many-spin model 


Single-spin model 


±10 


2.03-10-'^ K 


1.25-10-"' K 


±8 


4.20-10"'* K 


1.24-10-'^ K 


±6 


9.18-10"^ K 


1.20-10-'' K 


±4 


2.75-10"^ K 


2.27-10-^ K 


±2 


8.51-10-^ K 


7.01-10-" K 
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